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Classification of the universality exhibited by systems 
with macroscopic degrees of freedom, both at and away 
from equilibrium, is one of the main objectives that has 
been pursued in statistical physics ever since the ad- 
vent of scaling theory and renormalization group frame- 
work. The universality classes of nonequilibrium systems 
are far less understood, unlike those at equilibrium, in 
spite of having identified many genuine nonequilibrium 
classes like the absorbing phase transitions growing 
surfaces H, self-organized criticality 0], driven diffusive 
systems 4J, and so on. Constructing classes of infrared- 
stable field theories by taking a scaling limit of micro- 
scopic models is a formidable task even at equilibrium. 
Hence probing known field theories by various perturba- 
tions and following the induced instabilities, if any, is an 
alternative that can provide invaluable insights towards 
any classification. 

Near-equilibrium critical dynamics is extensively stud- 
ied and effectively captured by time-dependent Landau- 
Ginzburg (LG) models as categorized by Hohenberg and 
Halperin [5j. Recent studies have explored the effects 
of nonequilibrium perturbations on various dynamic uni- 
versality classes @, 0, B i, M, EE3- Th ey not only 
include perturbations about the LG energy functionals 
but also genuine nonequilibrium perturbations about the 
critical dynamics, ft is well established that the kinetic 
fsing systems of Model-A class (in Hohenberg-Halperin 
classification) are stable against local dynamic pertur- 
bations even if they violate detailed-balance condition 
provided the symmetries are preserved @, Il2j |. Bassler 
and Schmittmann (BS) further found that the spatially 
anisotropic perturbations in spite of not respecting the 
Zi symmetry cannot destabilise the dynamic class of non- 
conserved kinetic fsing models which are described by a 
single scalar field order parameter [7[. On the other hand, 
the detailed-balance violating perturbations turn out to 
be relevant in the conserved systems [H H, [Hj] . This natu- 
rally brings forth the issue whether the irrelevance of such 
perturbations pervades throughout model-A systems or 
is only restricted to its subset, like those describable by 



a scalar order parameter. To resolve this it is fruitful to 
understand the role of nonequilibrium perturbations on 
iV-component model-A systems. An incautious extension 
of the result of BS Q , namely the marginal irrelevance of 
the nonequilibrium anisotropic perturbations to model-A 
systems, may suggest that the A-component systems are 
also robust to such perturbations. Indeed such a miscon- 
ception seems to prevail in the community 0, [3] • 

The Gaussian fixed-point of equilibrium A-vector 
models in 4 — e dimensions is unstable to both Ising 
and 0(A)-symmetric perturbations and together they 
give rise to Ising, Heisenberg, and cubic fixed-points [l5( • 
The latter two are stable for N < N c and A > N c , re- 
spectively, where N c « 4 near 4-dimensions. In either 
of the case the Ising fixed-point is unstable. Now if we 
include nonequilibrium perturbations with model-A dy- 
namics then intuitively one may expect the dynamic class 
to be either model-A Heisenberg or cubic class or else 
some other A-vector nonequilibrium class depending on 
whether these perturbations are irrelevant or relevant. 
It is hard to imagine that the large-scale properties of 
the model could be governed by a kinetic Ising class. As 
we shall see, the nonequilibrium perturbations exert non- 
trivial effect on these models and make them cross over 
to a single component kinetic Ising class. We shall now 
describe how this crossover takes place. 

Consider the nonconservered A-vector model that is 
accorded with permutation symmetry while being sub- 
jected to spatially anisotropic dynamic perturbations, 
and described by the Langevin equation 



d t ip a {x,t) =J 7 a (tp(x,t))+r) a (x,t), 



(1) 



where 



T a ((p) = (V 2 - m 2 )ip a + 

+ J~] E abc d\\ ((pbVc) - J~] G abcd (p b ip c ip d , (2) 



be 



bed 



the indices a,b,c run over from 1 to A, and i] a (x,t) 
denotes the Gaussian noise with zero mean, and vari- 
ance (rj a {x,t)7] b (x',t')) ri = 2TS ab S{x - x')5(t - t'). The 
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bias couplings E a \, c distinguish one of the d dimensions, 
at ii, from its orthogonal space xj_, and introduce spatial 
anisotropy. Unlike the other terms in equation @ they 
spoil the symmetry of the dynamics (f> a i— ► — <f> a for any a . 
The permutation symmetry in this model is assured if the 
coupling constants are such that G-paPbVcPd = Gabcd and 
E-paVbVc = E a b c for all permutations V . Note that the 
above interaction terms are the most general marginal 
perturbations with permutation symmetry at d = 4. 
Terms like tp^dnipc can be set to zero without any loss 
of generality for this symmetry also demands E a b c to be 
symmetric in b and c. Renormalization will of course 
modify the above form of !F a {f>)- 

In the following, we essentially address the ques- 
tion: which of the interactions are consistent with a dy- 
namic N -vector field theory having a single characteris- 
tic length-scale in the large-time limit? To this end it is 
convenient to analyse in Martin-Siggia-Rose (MSR) for- 
16j . The MSR action of the Langevin equation 



the coupling constants do not satisfy the relation, 



malism 



(P) is given by 



TtpaVa 



(3) 



where, tp a refers to the auxiliary field, and the notation 
<Pa = <fa{ x i~t) and <f>a = tpa{, x i~t) is used. The form of 
T a as given in Eq. is not the most general action 
in concord with the permutation symmetry, as there are 
various other allowed terms: J2ab^P a Vb, J^abVa^b, and 
Sab ifiad^ipb- Some of these terms will appear when loop- 
corrections are taken into account in the effective action 
that describes the long-wavelength modes. Most impor- 
tantly it is the off-diagonal mass term J2 a b ^Pa^b, if gener- 
ated, that will introduce an extra length scale and turn 
some of the components of the iV-vector field noncriti- 
cal. Which of the interactions will avoid this scenario of 
dynamical breaking of permutation symmetry? 

Let us first examine the near-equilibrium pertur- 
bations {E a b c = 0) that are consistent with the 
./V-vector permutation symmetry. This symmetry 
in the dynamics will restrict the number of inde- 
pendent couplings to seven, which are denoted as 
{Gini, G1112, G1122, G1123, G1222, G1223, Gi234j-- The 
notation Gnu refers to those couplings G a bcd where all 
the indices b, c, d are same as a, and G1112 is used when 
one of the indices b, c, d is different from a, and so on. The 
one-loop correction to the quadratic part of the effective 
action due to these perturbations is given by 

5 G SW(&<p) =3AT f G -bcc ip a (x,t)tp b (x,t), (4) 

Jxt a,b,c 

where A is a regularization dependent constant, and 
in dimensional-regularization with minimal subtraction 
(MS) scheme A = m d - 2 T{-d/2 + 1)/(Att)^ 2 . Hence, if 



JY 



^ Gabcc — gSab 



(5) 



c=l 



for some arbitrary g, then the mass-matrix has off- 
diagonal elements. All these elements have the same 
value due to the permutation symmetry, and in turn 
render the matrix with two different eigenvalues, one of 
which is N — 1 degenerate. Therefore the symmetry gets 
dynamically broken at one-loop when the constraint is 
not satisfied. This is essentially a single constraint, and 
when explicitly expressed is given by 



G1112 + G1222 + (N — 2)Gi223 — 0, 



(6) 



for N > 2. For N = 1, the only nonzero coupling is 
Gnu. 

If the model further possesses a ip a — ► — ip a symme- 
try for each a, then there are only two allowed cou- 
pling constants, the cubic coupling Gnu and the O(N)- 
symmetric coupling G1122, and the constraint ((6]) is au- 
tomatically satisfied. This special case is the time- 
dependent Landau-Ginzburg model with O(N) symme- 
try, broken by cubic anisotropies. Equation JSJ) does not 
imply that G a b c d is symmetric in a and b. If G a bcd are 
chosen to be symmetric in all the indices, which is also 
a necessary condition for the system to reach equilib- 
rium, then there are only five independent G-couplings 
as G1222 = G1112 and G1 22.3 = G1123, and Eq. (J5]) reduces 
to the "trace condition" [151 ] . 

Now consider the bias interactions which are gen- 
uine nonequilibrium perturbations. The symmetry con- 
ditions, E-p a -pb-p c = E a bc for all permutations V 1 and 
E a bc is symmetric in the indices b and c, will allow only 
four E-couplings: {£111, -E112, -E122, ^123}- The one-loop 
correction to the MSR effective action due to the bias 
interactions is 



— BT I 2_j Ea, Cl c 2 E, 



(7) 

2b Lp a (x,t)d 2 tfb(x,t), 



a,b,ci,C2 



where B = (2 — d)m 2 A in MS scheme. Apart from the 
diagonal term, <p a d 2 <p a , the loop-corrections also bring in 

the term, J2b VadnVb, that is uncensored by the symme- 
try. We denote the ratios of these two terms with respec- 
tive to lfa^'j_ ( Pa term in the 1-loop effective action by p 
and p, respectively. For a generic choice of couplings the 
value of p will be nonzero, if not at 1-loop then at some 
higher order, unless the only nonvanishing couplings are 
Em and Gnu. Thus there are nonzero off-diagonal cor- 
rections to the propagator, (<pafb) 7^ for s/d, which at 
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two-loops leads to the following off-diagonal mass term, 
S eg SW=3At( - 1 =-4=) x 

^2 / G aia 2 a 3 ai'fia 1 (x i t)ip a2 {x,t). (8) 

„_ „_ „ . Jxt 



a\ ,d2 ,az ,a4 



At this order, the effective action also has a term propor- 
tional to SqS^ which is essentially a correction to g in 
equation ((5]). The off-diagonal corrections imply that 
in order to avoid a crossover from a dynamic TV-vector 
class the system has to satisfy yet another necessary con- 
straint, 



N N 



2 , / , Gabcd — g'Sab, 
c=l d=l 



(9) 



for some arbitrary g' . Using equation |5]) in the above 
expression leads to the explicit relation, 

2G 1122 +2(A-2) (G U23 + G 1223 ) + (7V-2)(^-3)G 1234 = 0. 

(10) 

Are there any consistent subsets of coupling constants 
satisfying the two constraints, ([6]) and (fT0|) , and closed 
under loop-corrections? Let us extend, for illustration, 
the earlier special case with only two G-couplings, Gnu 
and Gn 22 , which has the additional ip a — ► — tp a symme- 
try in the absence of bias. If we weaken this symmetry by 
compounding its action with xu — > —xu then the allowed 
E-couplings are Em and £122 • These E-couplings give 
a nonzero p, and hence the constraint (|10|) would require 
G\\2i = 0. But the given set of couplings is not closed 
under loop-corrections without Gn 22 interaction and is 
incompatible with a single length-scale theory. So the 
0(A)-symmetric perturbations make the critical dynam- 
ics of this model to eventually cross over to BS model. 
In general, in the absence of any other symmetry than 
permutation symmetry, the only consistent sets are the 
trivial one {Gnu, -Efn} an( i the set with all the cou- 
plings. The former is the BS model, while the latter is 
vulnerable to crossovers. 

Rephrasing the analysis, if we consider a permutation 
symmetric dynamic A- vector model, (right at the outset) 
with 

T a {<p) = (Vj_ + pd 2 - m 2 )Lp a + ^2pd 2 ip b 

b 

+ ^ E abcd\\ (<Pb<Pc) - ^ GabcdtPbtPctfid, (H) 



be 



bed. 



where the quartic couplings G a b c d violate any of the two 
constraints, ([6]) and (fTO)) . then the 1-loop corrections will 
make the system cross over from a dynamic N- vector class 
in 4 — e dimensions. For, when these constraints are not 
satisfied, one of the eigenvalues of the mass-matrix is dif- 
ferent from the rest. Therefore we can either have a one 



component critical dynamics, which is the BS model, or 
a dynamic (A — l)-vector model with bias. If we decide 
to tune the system to (A — 1)- vector model then it will 
again cross over when any of the two constraints, with A 
replaced by (A — 1), fail to hold. This result is likely to 
hold even in lower dimensions since the arguments lead- 
ing to the crossovers solely rely on the loop-corrections to 
a most relevant operator, namely, the off-diagonal mass- 
matrix. 

To be specific, if a bias is introduced in any system at 
0(A)-symmetric critical point in d < 4, even for N < 4, 
then we expect that it will cross over to the kinetic Ising 
class to which BS model belongs to. We have confirmed 
this prediction by the following numerical study. 

Consider a 0(2)-symmetric model on a 2-dimensional 
lattice described by the Hamiltonian 



n = - 



(n,m) 



E 



r + 2d 



(12) 



where n is the site index, (n, m) denotes sum over all 
nearest neighbor pairs, <j> n = {4> n ,ii<t>n,2) is a real 2- 
component vector field, and 4^ :— <f> n The dynamics 
of the field <j) n . a in the presence of a bias is governed by 
the following Langevin equation 



d 



on 

d(f> n ,a 



Edn 



,)+»M,a(t), (13) 



where T} n>a is the white noise with correlation 
{rin,a{t)tln',a'{t')) = 8 nn >6 aa >5(t - t'), and d\\ (0* )O ) := 
fin+i a ~ Yn-i ai where n + 1 and n—1 refer to the two 
nearest neighbors of n along a specified direction. In the 
unbiased case, E = 0, the steady state of Eq. (|T3")) is the 
equilibrium state with partition function 



-in 



(14) 



We have taken the two-dimensional square lattice to 
be of size L x L with periodic boundary conditions. The 
values of u and the bias E are set to unity, u = E = 1. 
Equation (| 13[) is then integrated numerically by employ- 
ing Euler method with At = 0.0025. The initial condi- 
tion is taken to be <j) n _ a — $a,i f° r an realizations. The 
system sizes of L = 2 6 , 2 7 , and 2 8 are considered, and 
the equilibration time is set to 20 000. After equilbration 
we measured the magnetization M. = ^2, n {4>n,i, 't'n^/L 2 
as well as M 2 := \M\ 2 and M4 := Mf at every 5 unit 
times, namely after every 2000 iterations with the above 
mentioned At, and then obtained the averages for all 
these quantities. The critical point r c is located using 
the Binder cumulant 



U L = 1 



(A/ 4 ) 
3<M 2 )i 



(15) 
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" trace conditions" their viable field theories are required 
to satisfy. The dynamic renormalization group analysis 
suggests, contrary to the widespread belief, that model-A 
systems can succumb to such perturbations. In general 
these perturbations are highly relevant and will make a 
dynamic TV-vector model whether in Heisenberg or in cu- 
bic universality class to cross over to kinetic Ising class. 
We confirmed this scenario by a numerical study on a 
two-component c/> 4 model in two dimensions. 
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FIG. 1: (color online) Finite size scaling collapse using data 
for L = 2 6 (square), 2 7 (triangle), and 2 s (circle). The up- 
per (lower) straight line whose slope is | ( — |) indicates the 
expected asymptotic behavior of the scaling function up to 
a multiplication factor. Inset: Binder cumulants as a func- 
tion of r for different system sizes as in the main figure. For 
comparison, the critical Binder cumulant for the Ising class is 
drawn as a straight line. 



The critical exponents, f3 and v, are found from finite 
size scaling by taking the scaling form for */ (M2) to be 

J(M^ = L-V v f{{r c -r)L l l v ). (16) 

The asymptotic behavior of the universal scaling function 
/ is given by 



as y 
as y - 



00, 
-co. 



(17) 



Our results are shown in Fig. [TJ The data collapse 
with the asymptotic behavior Eq. (fT7]) is in good agree- 
ment with the critical exponents of two dimensional Ising 



171 ] . The critical point is located 



model, fi = \ and v = 1 
at r c = -0.9545 ± 0.0010 as shown in the inset of Fig. Q] 
The value of the critical cumulant is consistent with that 
of the Ising model on a square lattice (~ 0.6107) [l8j |. 
Thus the model with dynamics Eq. ([T5]) clearly shows 
the order-disorder phase transition and exhibits critical 
behavior unlike its equilibrium counterpart which can not 
undergo such a transition in two dimensions f]~9j ] . The ki- 
netic Ising behaviour is actually what is expected from 
our analytic study. The fact that the model shows equi- 
librium behaviour further suggests that the bias is ir- 
relevant, and thus also reconfirms BS result [3] in two 
dimensions. 

In summary, we have studied the effect of spatially 
anisotropic nonequilibrium perturbations on noncon- 
served TV-vector models, and obtained the necessary 
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